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Approximation  of  Chronopotentiometric  Responses  by 
Orthogonal  Collocation 


B.  Stanley  Pons 

Present  Address: 
Department  of  Chemistry 
University  of  Alberta 
Edmonton,  Alberta  T6G  2G2 


Abstract 


High  accuracy  simulations  of  chronopotentiometric  responses  are 
demonstrated  for  several  electrochemical  mechanisms.  Derivation  of  the 
discretized  equations  for  simple  electron  transfer  and  the  EC  mechanisms 
are  given.  The  equations  for  other  mechanisms  are  tabulated. 


t 

I 

t 


The  advantages  in  using  the  orthogonal  collocation  technique  for 

the  simulation  of  second  order  partial  differential  equations  de- 

(1  2 ) 

scribing  electrochemical  diffusion  problems  have  been  reported  ' 

The  discretized  equations  have  been  derived  for  several  electro- 

(1-4 ) 

chemical  mechanisms  under  the  control  of  chronoamperometric 
and  cyclic  voltairanetric programs.  This  paper  presents  the 
discretized  equations  necessary  to  simulate  the  chronopotentiometric 
response  for  two  mechanisms.  Extension  of  the  technique  to  more 
complicated  mechanistic  schemes  is  shown  to  be  straightforward,  and 
the  appendix  contains  further  examples. 


Symbol 

[A];  [B];  etc. 
[A*] 


Ai,j'*  Bi,j 


A 

B 


cA;  cB;  etc. 

da;  db;  etc> 


E 

E° 

F 

i 

J 

K 

k 

n 

N 


13 


Definition 

Molar  concentrations  of  species  A,  B,  etc. 

Bulk  concentration  (initial)  of  species  A 

Discretization  coefficients^  ^ 

2 

Electrode  area,  cm 
Dt/Z2 

[A]/[A°] ;  [B]/[A°];  etc. 

Diffusion  coefficients  of  species  A,  B,  etc. 
*  °A  =  DB  “  DC'  etC* 

Kronecker  delta  (=  1  if  i  =  j;  =  0  if  i^j) 
Potential 

Standard  potential 
Faraday 

Current,  amperes 
Elements  of  the  Jacobian  matrix 
Homogeneous  equilibrium  constant 
Homogeneous  rate  constant. 

Number  of  electrons  transferred 
Order  cf  polynomial  used 


Symbol 


Definition 


R 

T* 

T 

T 

t 

X 

X 

Xi 

Z 


Gas  Constant 
Temperature 
Time ,  seconds 

Chronopotentiometric  transition  time 
=  T/t 

Perpendicular  distance  from  electrode  surface  (cm) 
X/Z 

Roots  of  orthogonal  polynomial  (Legendre) 

Distance  in  solution  from  electrode  surface  where 
diffusional  effects  are  negligible  during  time  of 
experiment 


Simple  Reversible  Charge  Transfer 

ne~ 

A  ^  B  <  1) 

This  mechanism  is  considered  fof  the  planar  diffusion  case 
where  initially  only  species  A  is  present.  At  t  =  0,  a  constant 
current  i  is  applied  to  the  cell,  and  the  potential-time  plot  is 
recorded,  noting  the  transition  time  t.  The  response  is  given 
analytically  by^ 


E  =  E° 


RT* 

nF 


*5 


(  2) 


Although  this  solution  is  fully  adequate  for  describing  the 
response  after  corrections  for  charging,  derivation  of  the  discretized 
collocation  equations  will  be  given  so  as  to  exemplify  the  procedure 
that  is  followed  for  developing  the  approximation  of  the  responses 
of  more  complicated  mechanisms. 


The  diffusion  equations  and  boundary  conditions  describing  this 
experiment  are 


SIA)  _  _  S2(A] 
2 

6T  A 

6  IB]  _  62IB] 

«T  B«X2 


(  3) 

(  4) 


tAIx,0  '  l«-,t  ”  I*’! 

(  5) 

>®>x,0  *  *  0 

(  6) 

a* 

11 

a 

(  7) 

D  -  db 

(  8) 

/«IA]\  -  (6  [B]\ 

\  6X  /X*=  0  B^  5X  yx=o 

(  9) 

i  -  nFADA^  6X  jx=0 

(10) 

£g 

(11) 

Insertion  of  the  dimensionless  variables  for  concentration, 
distance,  and  time  into  (3)  and  (4)  gives 


6c  a  D  [A0  ]  u  CA 


r  ft 


226x2 


(12) 


|A°)k’B  D  [  A0  ]  6  Cfi 

"  7I?- 


T  ft 


(13) 


and  after  simplifying. 


0 


et  t  x2 


(14) 


(15) 


The  boundary  conditions  (5) -(7)  are  treated  similarly  and  become: 


cA(x,0) 

=  1 

cA(l,t)  =  1 

(16) 

°B  (x,  0) 

*  0 

cfi(l,t)  «  0 

(17) 

ft 

X 

II 

o 

II 

Ji !b\ 

/x1=0 

(18) 

The  discretization  equations  for  first  and  second  order  ordinary 
differentials  have  been  given  previously ^ ,  and  are: 


dy (x,t) 
dx 


N+2 

1  A.  .  y  (x  .  ,t) 
j=l  K'3  3 


(19) 


d^y  (x,  t) 


dxd 


xi 


N+2 

1  B.  .  y  (x  . ,  t) 
j=l  1,3  3 


(20) 


when  evaluated  at  the  roots  of  an  orthogonal  polynomial.  (See 
computational  section.) 


Application  of  (20)  to  (14)  and  (15)  at  the  roots  of  a  Legendre 
polynomial  yield 


3t“  x±  Bj!1Bi,j 


U^B  _  6  ;  Bi  i  cB(xVt} 
dt^  x.  ^  1/3  B  3 


Expanding  (21)  and  (22)  partially,  we  have 


dt  x.  e*Bi,lCA(0/t)  +  Bi,N+2CA(lft)  +  ,£,Bi, jCA (xj rt) 1  (23) 

1 


5tT  x.  “  BlBiacB(0,t)  +  B.#N+2cB(l,t)  +  I  Bi#jcB(*jft))  (24) 

,  ™  J 

Insertion  of  the  boundary  conditions  (16)  and  (17)  yields 


®Txt  ■  »[Bi,lCA<0't)  +  Bi ,N+2  +  j£lBi,jCA(l'yt>1 


at5**  *  +  jf1Bi,jcB(xrt)1 


We  need  an  expression  for  cA(0,t)  and  cB(o,t)  explicitly.  From  (20) 
applied  to  the  boundary  flux  condition  (18) ,  we  have 


53T  Xj-0  ’  .?lAl 


35“  *ro"  ,*,*), 1,cn (xi't’ 


Now,  equation  (10)  is  also  made  dimensionless  as  follows. 


Combining  this  expression  with  (27)  and  (28),  we  have 


N+l 

c *  ( 0 , t )  =  Q  -  I  a.ca(x.  ,t) 
A  j=2  3  A  3 

N+l 

Ct>  (0,  t)  =  R  -  I  a  ,c_  (x . ,  t) 
B  j=2  3  B  3 


where 


Q  = 


*[ 


it** 


Al,l  |nFAlA°jD*B* 


A1 ,  N+2J 


R  =  _I_  _  ±lL  1 

Al,l  LnFA[AelD,J3,5  J 


,  and 


-  & 
3  1,1 


Since: 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


Substitution  of  (31)  and  (32)  back  into  (25)  and  (26),  we  have 


9 


where 


(37) 


dcB 

dt 


Ti  +  6 


N+l 

I 

j=2 


,t) 


Si  =  PlBi,lQ  +  Bi,N+ 21 


(38) 


T. 

l 


6  B1(1R 


and 


(39) 


+ 


(40) 


Simultaneous  solution  of  the  2N  differential  equations  by  the 
integration  subroutines  described  previously  give  accurate  and  fast 
approximations  to  the  consentration  profiles  of  A  and  B  as  a  function 
of  time. 

The  desired  response,  potential  as  a  function  of  time,  is  given 
then  by  equation  (11)  in  the  form 


with  cfi(0,t)  and  cA(0,t)  being  supplied  at  each  time  integration 

step  by  equations  (31)  and  (32).  A  comparison  of  some  typical 

E-t  points  from  equation  (2)  and  from  equation  (41)  is  given  in 

Table  1,  and  the  simulated  plot  shown  in  Figure  1.  Only  a  sixth 

order  Legendre  polynomial  (N  -  6)  was  used.  Total  computational 

time  for  1000  E-t  points  was  0.862  seconds  (see  Computational  sec- 

.  .  (7) 

tion)  .  8  may  either  be  determined  explicitly  ,  or  fixed  arbitrani 

(2 ) 

and  integrated  by  splines  '. 

EC  „  Mechanism 
rev _ 

The  extension  of  the  matin  m.j tics  to  include  kinetic  react  ions 


is  immediate.  Consider  (under  the  same  experimental  conditions)  the 
ECrev  mechanism: 


k  , 

with  K  =  ,  and 

K1 


(43) 


*  -  (k^  +  k_1)x 


(44) 


The  diffusion-kinetic  equations  and  boundary  conditions  de¬ 
scribing  the  system  are: 


6  [A]  6*  [A] 

ST  A  {X2 

(45) 

+  k_x  £CJ 

(46) 

-  k_x  tci 

(47) 

IA'x,o  “  ^a^«,t  “  0 

(A*) 

(48) 

[B1X,0  -  fBl-.T  -  0 

(49) 

IC1X,0  "  IC>«,T  '  0 

(50) 

n  _  _D  /d  [B]\ 

DA\^>ryx=o  db\i"it7x=o 


(51) 


.1 


dCg(t) 


dt 


N+l 


x.  =  S[BirlcB(0,t)  + 


t)]  - 


Ml+K)”1!  CgCxj^/t)  -  Kcc(xi,t)] 


(62) 


dcc(t) 


dt 


N+l 

=  6  IB.  ,cr  (0,t)  +  I  B.  .cr(x.,t)]  + 
x£  +»+  j_2  +  »3  ^  3 


A  (1+K) -1  [cB  (x^, t )  -Kc^x^t)] 


(63) 


The  flux  relation  (59)  is  identical  to  the  simple  reversible  charge 
transfer  case,  and  is  discretized  in  precisely  the  same  manner  to 
yield  equations  (31)  and  (32) .  The  flux  relation  for  species  c  is 
equation  (60) ,  and  is  discretized  as 


(—) 

\6x  /x=0 


N+2 

=  0  =  E  A,  -c  (x.,t) 
j=l  c  3 


(64) 


N+l 

A1,1CC (0,t*  +  Al,N+2GC^1,t^  +  Al, jCC*Xj' V 

J”2 


N+l 


=  A.  .0^(0,^  +  E  A,  .c_(x.,t) 
*»A  j_2  A*J  ^  J 


(66) 


or 


N+l 

c_(0,t)  =  E  a.c  (x.,t) 
c  j=2  3  C  3 


(67) 


Substitution  of  (31),  (32)  and  (67)  into  (61),  (62),  and  (63)  gives 


dcA 
_ A 

dt 


N+l 

+  B  E  b. 
j"2  J' 


jCA(Xj 


t) 


(68) 


i pr 


dc 

at 


B 


N+l 


-1 


x  -  Ti  +  si  bi,jcB(xj,t.)  "  X<l+K)  I  cB(xi,t)-Kcc(x.,t)]  (69) 
*  j~2 


ffc 

dt 


N+l 


-1, 


=  3  E  b.  c  (x.,t)  +  X  (1+K)  [c_  (x.  ,  t)  -Kc-,  (x . , t) ] 

j_2  ^/J  J  o  i  -  v.  i 


(70) 


Again,  simultaneous  solution  of  equations  (68)- (70)  provide 
the  time  dependent  concentration  profiles.  The  chronopotentiometric 
response  for  this  mechanism  is  still  given  by  equation  (41)  with 
the  term  in  parentheses  being  supplied  by  equations  (31)  and  (32). 
The  concentration  terms  in  (31)  and  (32)  however,  are  now  calculated 
numerically  from  equations  (68)-(70). 

The  results  of  the  simulation  of  this  mechanism  was  compared 
to  the  response  given  by  the  analytic  solution ^ 


E  =  E° 


RT' .  T 
"F  ln  ti>_t 


,  rt*.  r 

h  +  nF  ln' 


Ll+k 


m»jK  erf  [  (k1+k_1 )  ^T*5 ]-j 

2  (1+K)  (k-+k  ,)V  J 


(71) 


The  comparison  of  the  two  solutions  shows  that  the  simulated 
results  are  easily  maintained  within  a  relative  error  of  0.01% 
when  3  is  chosen  by  the  method  described  by  Rieker  and  Speiser^. 

Simulations  for  more  complicated  mechanisms  are  formulated 
simply  by  the  same  procedures  used  to  arrive  at  equations  (68)  —  (70) . 
Under  similar  boundary  conditions,  the  results  are  quite  easily 
obtained . 

Computational 

All  simulations  were  performed  on  the  Amdahl  470V7  computer  at 


the  University  of  Alberta.  The  main  integration  subroutine  for  the 
system  of  coupled  first  order  ordinary  differential  equations 


((36) -(37)  or  (68) -(70))  was  the  same  third  order  semi-implicit 

(2) 

Runge-Kutta  method  as  described  previously  for  "stiff"  equations. 

The  subroutine  is  very  general  in  usage,  and  calls  three  external 

(8 ) 

subroutines  furnished  by  the  user  :  the  first  supplies  the  algorithm 
for  the  right  hand  side  of  the  time  derivatives  of  each  concen¬ 
tration  term  (equations  (36)-(37)  or  (68)-(70)),  the  second  is  a. 
simple  output  routine  for  displaying  the  results  and  is  called 
automatically,  and  the  third  is  a  user  supplied  matrix  of  the 
Jacobian  of  equations  (36)  —  (37 )  or  (68)-(69),  necessary  during  the 
integration . 

The  elements  of  the  Jacobian  matrix  are  simply  the  derivatives 
of  equations  (36) -(37)  or  (68) -(70)  with  respect  to  each  x^.  The 
matrix  elements  of  the  Jacobian  for  the  simple  reversible  electron 
transfer  mechanism  then  are  given  by: 


and 


for  species  A 


for  species  B 


These  are  given  explicitly  by: 


for  species  A,  and 


+ 


N+l 
Z  b 
j“2 


i,  jCA 


J 


i# 


+ 


N+l 
Z  b 

3  =  2 


i ,  jCB 


(72) 


(73) 


(74) 


(75) 


for  species  B,  or  finally, 


(76) 


Pi,j 


"Bi,lAl,j 

*1.1 


+  B 


for  both  species. 

Similarily,  the  matrix  elements  of  the  Jacobian  for  the  ECrev 
mechanism  are  given  by: 


J,  .  •  equation  (76) 
if  J 


for  species  A,  and 


(77) 


— B .  A  » 

J,  ^  =  — 1^1  1,3  +  B,  4  -  XK(l+K)"1cB(xi,t)6 _ 


If3 


Mfl 


irD 


13 


(78) 


for  species  B,  where  6^  is  the  Kronecker  delta,  and 


—B.  .A.  . 

J.  .  »  — +  B  .  -  XK(1+K)”1c_(x.  ,t)6  .  . 
if  3  A11  if 3  Ci.  13 


for  species  C. 


(79) 
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Table  1 


E/volts,  E/volts 


Time/s 

simulated 

exact 

0.0023 

0.0771 

0.0771 

0.0107 

0.0556 

0.0556 

0.0418 

0.0351 

0.0351 

0.1188 

0.0167 

0.0167 

0.3358 

-0.0048 

-0.0048 

0.5055 

-0.0227 

-0.0227 

0.9100 

-0.0752 

-0.0752 

0.9940 

-0.1241 

-0.1240 

Comparison  of  simulated  results  with  exact 
(equation  2)  for  simple  reversible  charge 
transfer,  i  =  0.269  ma,  t  =  1.010  s,  D  = 
10“5  cm2/sec,  [A°J  «  1  x  10“®  mole/cm2, 

T  *  25°C,  A  =  1.000  cm2#  linear  diffusion, 
quiet  solution,  E°  =  0.000  V. 


Calcomp  digital  plot  of  simulated  chronopotentiogram 
for  simple  reversible  charge  transfer  mechanism.  See 
Table  1  for  experimental  parameters. 


Listed  below  are  the  discretized  equations  for  several  other 


commonly  occurring  electrochemical  mechanisms.  As  above,  it  is 
assumed  that  the  diffusion  is  to  a  planar  electrode  in  quiet 
solution,  and  species  A  is  the  only  electroactive  species  present 
initially. 

Catalytic  Mechanism 

ne~  k 

A  7—**  B  - >  A  +  C  (80) 


dCj 

dt 


lx.  =  Si  + 
1 


N+l 

8  E  b.  .c  (x.,t)  -  kxcB(x.,t) 
j_2  °  3  a  D 


(81) 


dcB 

dt 


N+l 

T.  +  8  E  b.  .c  (x. ,t)  -  kxcB (x . , t) 


j=2 


(82) 


For  the  E-t  profile,  equation  (11)  is  used  along  with  equations 
(31)  and  (32) ,  the  unknown  concentration  terms  being  furnished  by 
the  simultaneous  solution  of  the  2N  equations  (81)  and  (82) . 


Dimerization 


A 


ne 


B 


We  let 


X  «  (lA*Jk.  +  k_^)x,  and 


(83) 


(84) 


K  « 


(85) 


Then 


3T—  »  S.  +  0  Z  b.  .c  (x.  ,t) 
at  x  ^2  +  »J  A  J 


dcB  N+1  -1  2 

at"  =  Ti  +  P  E  b± #jCB(Xj,t)  -  X  (1+K)  ■L(cB(xi,tr  -  Kcc(xi,t)) 
j=2 


fic  N+1 

St2"  642bi,jCC(Xj,t>  +  X(l+K)-1(l,cB(xi,t)2  -  Kcc(xi,t)) 


Once  more,  equations  (11) ,  (31) ,  and  (32)  are  used  for  the  E-t 
profiles.  Note  that  in  this  case,  3N  equations  must  be  solved 
because  cB(0,t)  depends  on  the  concentration  of  the  C  species. 


Second  Order  Reaction 


ne 

a  r=±  b 


B  +  A  rz±  c 
k  , 


X  and  K  are  defined  as  in  equations  (84)  and  (85) . 


-  Si  +  p  l  jCA(x..  ,t)  -  X  (l+K)“1(cA(xi,t)cB(xi,t)  - 

j»=2 


Kcc(xi,t)) 


dc  N+1  , 

dT  Ti  +  P^2bi,  jCB<xj  ,t}  "  *^1+K)  (cA(xif  t)cfi(xi,t) 


Kcc(xj ,t) ) 


